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NOMENCLATURE 


flow coefficient, nozzle mass flow/1-D mass flow 
specific Impulse 

specific impulse for an isentropic equilibrium 
expansion 

Mach number 

radial coordinate. r-o at axis 

throat radius 

throat radlus/r* 

downstream throat radius/r* 

upstream throat radius/r* 

axial coordinate. x«o at throat plane 


Greek 

e 

Y 

e 

’^DER 

^TDK.BLM 

e 


angle between shock and gas streamline, 

shock strength (see page 4) 

ratio of specific heats 

nozzle expansion ratio 

Distributed Energy Release efficiency 

aerodynamic efficiency, IncTUdlng chemical 
kinetics 

angle between gas streamline and axis 
angle between surface and axis 



^sp 

^sp 

^theoretical 

M 

r 

» 

r 

R 

«d 

R 

u 

X 
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ABBREVIATIONS 


BLM 

DER 

JANNAF 

LRC 

MCC 

M0C 

0DE 

0DF 

0DK 

RRC 

SEA 

SSME 

TDK 

VNAP2 


Boundary Layer Module, computer program 

Distributed Energy Release 

Joint Army-Navy-NASA-Air Force 

Left Running Characteristic 

Main Combustion Chamber 

Method of Characteristics 

One-Dimensional Equilibrium Expansion 

One-Dimensional Frozen Expansion 

One-Dimensional Kinetics Expansion 

Right Running Characeristic 

Software and Engineering Associates, Inc. 

Space Shuttle Main Engine 

Two-Dimensional Kinetics, JANNAF computer program 
Viscous Nozzle Analysis Program, Version 2 
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1 . INTRODUCTION 


—In -this study the Two-Dimensional Kinetics (TDK) computer 

program that has recently been developed for the NASA under 

contract NASS- 35931. was used to predict the performance of a 

"Large Throat Configuration" of the Space Shuttle Main Engine 

(SSME)^^\ Calculations indicate that the current design SSME 

contains a shock wave that is induced by the nozzle wall 
( 2 > 

shape . In the "Large Throat" design an.even stronger shock 

wave is predicted. Because of the presence of this shock wave, 

earlier performance predictions that have neglected shock wave 

effects have been questioned. In this study the JANNAF thrust 

chamber performance prediction procedures given in Reference 3 

have been applied. The analysis includes the effects of two 

dimensional reacting flow with a shock wave. The effects of the 

boundary layer with a regenatively cooled wall are also 

included. This analysis has been made possible by the 

development of a new shock wave version of the TDK computer 
(1 ) 

program 

In order to substantiate the calculations, a computer 
program' developed by Prof. Hoffman at Purdue University was also 

( H ) 

used . The results were run independently at Purdue with 

input data supplied by SEA, Inc. The Purdue computer program 
can compute axially symmetric supersonic nozzle flows with an 
induced shock, but is restricted to flows with a constant ratio 
of specific heats, Y . Thus, the TDK program was also run with 
this assumption and the results of the two programs were 
compared . 
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Nozzle Geometry and Operating Conditions 


The Space Shuttle Main Engine (SSME) utilizes a thrust 
chamber that is constructed in two sections, 1) a Main 
Combustion Chamber (MCC) assembly, and 2) an expansion nozzle 
assembly. The MCC current design is shown in Figure 1 . The 
nozzle assembly is bolted to the MCC at a nozzle attach flange. 
The nozzle expansion ratio at the attachment position is 5:1. 

The "Large Throat" version of the SSME is constructed by 
modification of the MCC. The nozzle assembly is not modified. 
The current design and "Large Throat" design MCC geometries are 
compared in Figure 2. It can be seen that the enlargement of 
the nozzle throat reduces the nozzle expansion ratio from 77.5 
to 69.5:1. The chamber pressure is reduced from 3285 (109$ 
power level) to 3010 psia. The distance from the nozzle throat 
to the nozzle exit is very slightly decreased. There is a 
considerable difference in throat radius of curvature, both the 
upstream and downstream values, between the two designs. 


Nozzle geometry and operating conditions^ for the SSME 
current design nozzle (109$ power level) were obtained from 
NASA/MSFC and are given in Figure 3, which is an annotated 
listing of the input data. Definitions for all of the input 
data items are given in Reference 1 . Those input items which 
were changed in order to run the "large throat" SSME are listed 
in Figure 4. The wall temperature profile for these engines is 
shown in Figure 5. 


The calculations were made assuming that 1) there is perfect 
mixing so that the overall engine mixture ratio of 6.05485 is 
achieved, and 2) the system enthalpy is determined from tank 
conditions. These assumptions are equivalent to assuming that 


the distributed energy release efficiency, n 


DER 


is unity. 
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Results were then calibrated using measured data. That is, an 
actual value of n^g^was estimated from 



The performance predictions results are presented in the 
next section. 
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PERFORMANCE PREDICTION RESULTS 


Performance prediction results are summarized in Figures 6 
and 7 for the current design and "Large Throat" designs, 
respectively. (See Appendix A for an explanation of these 
summary sheets). The measured Igp for the current design is 
known to be ^52.6±-1 seconds. The distributed energy release 
efficiency for the engine is found from equation (1) as follows: 

I - 458.0 sec. 

TDK.BLM 

"TDK, BLN - ■'58.0/465.1 - .9847 

"DEB - = 

Thus, it is predicted that the distributed energy release 
loss for this engine is 5.4 seconds. This value is consistent 
with the value of 5.3 seconds predicted by R. Carroll^^^ using a 
full distributed energy release analysis. 

As can be seen from Figure 6, the maximum shock strength, 6, 
predicted for the current design nozzle is .29. This value 

is taken at the intersection of the last Left Running 
Characteristic (LRC) surface and the shock surface. Downstream 
of this intersection the shock can not effect thrust chamber 
performance. Shock strength is defined here as 

e - (Pp - ) / P, 

where 

P.| is the pressure in front of the shock, and 

Pp is the pressure behind the shock. 
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For the large throat design the maximum shock strength is 

found from Figure 7 to be 1.13. The difference in theoretical 

performance between the two engines is (ISP (ODE) values from 
Figure 6 and 7) 

465.13 - 463.26 - 1.87 seconds. 

The difference in predicted performance is (ISP (TC) *nDER from 

Figure 6 and 7) 

452.6 - 449.1 - 3.5 seconds, (with nDER = .9882) 

which is larger than expected. From the analysis presented here 
it appears that the lowered performance predicted for the "Large 
Throat" SSME design is almost entirely due to the effect of the 
shock wave calculation. Factors that influence the performance 
calculations are discussed in the next section. 
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3. FACTORS EFFECTING THE PERFORMANCE PREDICTION RESULTS , 


An examination of the computed results indicate that the 
lowered performance predicted for the "Large Throat” SSME is a 
result of the shock wave calculation. In order to minimize the 
numerical error contained in the calculation, an extremely fine 
mesh spacing was used. The number of points in the initial data 
line was set at 250. Each Method of Characteristics (MOC) flow 
field calculation was carried out using approximately 100,000 
points . 

Two phenomena were examined in detail; 1) the interaction 
of the shock wave with the chemical kinetics of the flow, and 2) 
the transonic flow field. Each is discussed below. 

Shock Wave and Chemistry Interaction 

Across the shock front, the translational degrees of freedom 

* 

of a gas are equilibrated within a few molecular collisions . 
The chemical composition of the gas is frozen during this time 
since thousands of molecular collisions are required for even 
the fastest chemical reactions. Thus a chemical excitation 
process can occur immediately behind the shock front. Since the 
kinetic rate equations are notoriously stiff, it is necessary 
that a fully implicit integration method be used (see Ref. 1, 
Section 3.1.1). 


* The present analysis assumed that both the translational and 
vib.—rot. degrees of freedom are equilibrated across the 
shock . 
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In order to test the accuracy of the TDK integration 
procedure immediately behind a shock point, a streamline was 
selected froitt the "Large Throat" SSME TDK analysis. A pressure 
defined ODK calculation was then carried out for this streamline 
starting at the back side of the shock. The shock strength was 
1.05, so that there was a considerable change in pressure and 
temperature across the shock. 

A comparison of the TDK and ODK streamline calculations are 
shown in Figure 8. Each TDK mesh point is shown. The ODK 
integration step is approximately 1000 times smaller than the 
TDK mesh size. Never the less, the two results cannot be 
separated on the scale shown in the figure. Velocity vs path 
length was chosen since velocity is a direct measure of specific 
impulse (for one-dimensional flow, optimum Igp - V/g). A 
computer generated plot repeating this result, but with higher 
resolution, is given in Figure 9. It can be concluded that for 
this case, TDK is giving an accurate result. 

A one-dimensional frozen (ODF) calculation was also carried 
out for this streamline. There were no differences between the 
ODF and- ODK results, i.e., the frozen and equilibrium paths were 
found to be the same. Thus, in this flow regime there can be no 
dissociation of HgO or due to a shock of this strength. It 
follows that no significant energy can be bound up in 
dissociated species causing a performance loss. The shock wave 
and chemistry interaction is founts to be negligible. 

The Transonic Flow Field 


The most important parameter effecting the transonic flow 
field is the radius of curvature immediately upstream of the 
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nozzle throat. This value was reduced from • 1 in the 

current design to » .494 in the "Large Throat" design (see 

Figure 2). The modified Sauer transonic analysis used by the 

TDK program is based on perturbation theory (see Reference 1, 

Section 2.4) and a value of R^ - .494 is at the lower limit for 

which the analysis has been validated. Because of this fact, 

the analysis was re-examined in order to assess its accuracy for 

calculating the MOC start line for the "Large Throat" SSME. In 

particular, results obtained from the analysis were compared 

( 7 ) 

with results obtained from the VNAP2 computer program . Two 

cases were examined in some detail. 

The first case examined was the nozzle studied by L. Back, 
Reference 8, for which ample experimental data is available. 
The Back case is for air flow through a nozzle with a 
small, R^ - .625, radius of curvature throat. Results of the 
analysis -are shown in Figure 10. It can be seen that both 
methods are in excellent agreement with the data in the region 
near the wall. Near the flow axis, the TDK analysis is in 
closer agreement with the data than the VNAP2 analysis. In 
general, VNAP2 predicts flow accelerations along the nozzle axis 
that are greater than seen experimentally (see Reference 8 for 
further data). 


The sedond case examined is the "Large Throat" SSME design. 
Results of the VNAP2 analysis are shown in Figure 11. Included 
in Figure 11 is the isomach attached to the throat minimum point 
as calculated for this case by the TDK modified Sauer transonic 
analysis. It can be seen that the agreement between this 
isomach and the corresponding isomach computed by VNAP2 is very 
similar to that measured by Back, et.al.. Figure 10. A grid 
spacing of 101x34 was used for the VNAP2 runs. The VNAP2 
solution grid for the "Large Throat" SSME case is shown in 
Figure 12. 
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For the prediction of rocket engine performance it is 
important to be able to accurately estimate the steady state 
mass flow rate. This parameter is both difficult to measure 
experimentally or to predict theoretically for nozzles with 
small throat radius of curvature. The mass flow coefficient, 
Cp , can be shown to be rather sensitive to the flow angle 
distribution within the nozzle. Consider the isomach attached 
to the throat minimum point for the Back case shown 
in Figure 10. For this suface (M - 1.3*»6) the depends only 
on the flow angle distribution, 0(r). The will be maximized 
if the flow is normal to the surface (although this violates the 
wall boundary condition). On the other hand if the flow angle 
is parallel to the nozzle centerline, then the Cp is very much 
reduced. Figure 13 gives computer values for Cp for these two 

cases (e » (}) - ■! and 0-0). Also given is the Cp corresponding 
to the flow angle distribution computed by the TDK modified 
Sauer method. This computed distribution has then been 
perturbed + 10$ and values of Cp determined. It can be seen 
from Figure 13 that a 10$ perturbation in flow angle 
distribution gives significant differences (.5 to .8$) in 
computed mass flux. Results obtained for Cp and for the isomach 
attached to the nozzle throat minimum point are shown in Figure 
14 . 

From the results discussed above it is concluded that the 
TDK modified Sauer analysis is adequate for the present 
application. 
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4. 


COMPARISON OF SHOCK RESULTS WITH OTHER COMPUTER PROGRAMS 


The current design and "Large Throat" design nozzles were 
analyzed using a method of characteristics computer program with 

( ii ) 

shock tracing that was developed by Hoffman . This MOC 
computer code can only analyze nozzle flows where the ratio of 
specific heats, Y, is constant . A value of Y-1.24 was selected 
for the analysis since this closely approximates the equilibrium 
value found in the region of the shock, i.e., at T ■ 2500 °R. 
These same two cases were also run using the constant Y option 
of the TDK computer program. A comparison of results obtained 
are presented in Figures 15 through 20. 

The nozzle contour for the current configuration of the SSME 

is shown in Figure 15. An induced shock wave is found by both 

(4) 

the TDK computer program and the Hoffman computer program , 
which is indicated in the Figures by the initials JDH. The 
corresponding plot for the "Large Throat" configuration of the 
SSME is shown in Figure 16. 

A plot of shock strength along the shock path for the 
current configuration of the SSME is shown in Figure 17. The 
corresponding plot for the "Large Throat" configuration of the 
SSME is shown in Figure 18. 

From the above figures it can be seen that the shock begins 
as a weak disturbance. The shock path follows a Right Running 
Characteristic (RRC) surface until x/r* - 4 for the current 
design nozzle, and x/r* - 3 for the "Large Throat" design 

nozzle. The shock wave then gains in strength until the last 
LRC surface in the nozzle is reached. Downstream of this 
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position the shock cannot effect nozzle performance. 

It was found that the shock strength calculation is 

sensitive to the input parameters. A sensitivity study was 

therefore carried out by J. D. Hoffman using the computer 

program described in Reference 4. The number of points on the 

initial data line was varied by doubling the spacing as follows: 

21, 41, 8l , and I6l. The computed values of shock strength vs 

position are shown in Figure 19. The Hoffman program uses the 

( 9 ) 

modified Hall method for transonic analysis in constructing 

the initial data line. In order to determine the effect on 

shock strength of the transonic analysis, an initial data line 
calculated by TDK using the modified Sauer method was input 
into Hoffman's program. The results are also shown in Figure 19. 

A sensitivity’ study was also carried by SEA using the TDK 
computer program. The purpose of the study was to examine the 
sensitivity of the performance .calculation to the pressure 
calculation behind the shock front. It was found that a 1 Jl 

error in pressure corresponds to roughly .75 seconds in specific 
impulse . 

The Lockheed HOC program^^^^ was also Used to calculate the 
induced shock wave in these nozzles. Both Lockheed Huntsville, 
and Rocketdyne, used their version of the program to perform 
the calculations. Equilibrium chemistry was assumed. The 
results of all of these calculations are summarized in Figure 
20 . 
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5. 


CONCLUSIONS 


As a result of the studies described in this report, the 
following conclusions are given regarding the performance of the 
SSME engine designs. 

5 . 1 Chemical Kinetics 

The effects of chemical kinetics are not important. The 
difference in predicted specific impulse between the ODE and ODK 
analysis is only 2/10 seconds. It is also shown that the shock 
waves found in the nozzles do not produce a kinetic loss. 

There is no significant difference between equilibrium and 
frozen gas properties in the region of the shock. A Mach shock 
would be required to activate the flow reaction chemistry. 

5 . 2 Transonic Analysis 

Three different transonic methods were used to obtain start 

lines for the MOC calculations, 1) modified Sauer^^\ 2) 
( 9 ) 

modified Hall , and 3) a time dependant finite difference 
(7) 

method . Each method produced a somewhat different shock wave 
solution. It is not possible to determine which transonic 

analysis is the more accurate. The modified Sauer analysis 
gives the best prediction of measured data'’ . The VNAP2 
analysis, however, is based on more advanced methods. A nozzle 
design should be suspect if it shows a significant shock wave 
based on calculations using any one of these transonic methods. 
If the shock wave becomes sufficiently strong, e.g., B > 1, then 
according to the TDK computer program the nozzle performance is 
significantly degraded. 
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5.3 Shock Wave Prediction 

As shown in Figure 20, nozzle flow field calculations were carried 
out independently by several contractors. Each analysis predicted the 
existence of a shock wave in both the current and "Large Throat" SSME 
nozzle configurations. The SEA prediction gave a shock strength for 
the "Large Throat" nozzle that was 2.7 times as strong as in the current 
SSME design (see column 1 of Figures 6 and 7). Changes either in the 
transonic method, or in the MOC grid spacing were found to affect the 
shock strength as much as 50%, however, the predicted performance of 
the nozzles was little affected. 

5.4 Thrust Chamber Performance Prediction 

The analysis implies that for the current design SSME thrust 
chamber, the performance loss due to Distributed Energy Release is 
more than 5 seconds of specific impulse. If this is true, then 
improvement of the injector design could recover most of this 
performance loss. 

The TDK computer program predicts a substantial* performance 
loss for the "Large Throat" SSME configuration, as compared to the 
current design. The loss appears to be primarily due to flow divergence, 
rather than shock wave effects. Because of the small downstream throat 
radius (Rj = .2) of the "Large Throat" SSME a finely spaced MOC grid is 
required near the wall to obtain an accurate performance prediction. The 
numerical accuracy of the performance predictions require further 
investigation. 

The calculations indicate that a weak. shock has little effect on 
performance, but a stronger shock (3 1) can have an important effect 

on performance. The "Large Throat" SSME configuration appears to be 
near the threshold where shock wave effects become important. 


*. Loss due to nozzle aerodynamics = 3.5 - 1.9 = 1.6 seconds. 

See Section 2. (Aerodynamic Loss Includes: Divergence with Shocks and 

Boundary Layer Effects) 
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6. RECOMMENDATIONS 


It is recommended that the SSMET MCC be recontoured to reduce 

or eliminate the shock wave predicted for the "Large Throat" 

SSME design. This has recently been done by Rocketdyne 

by replacing the contour downstream of the nozzle throat with a 

( 13 ) 

shape designed with the aid of a Rao method computer 

program. The revised contour was run by SEA, Inc. using the TDK 
computer program, which showed the nozzle to be shock free. 
The predicted nozzle performance was increased by .3 seconds. 


With regard to a projected "Large Throat" SSME test program, the measurement 
of the shock location in the exit plane is recommended. A comparison of the 
predicted with the measured value will contribute significantly to verify 
the current analytical calculation method. 
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TITLE SSME - COOLED WALL , CURRENT DESIGN, 109X 
DATA 
SOATA 

OOE * 1, OOK * 1, TDK » 1, 8LM » 1, IRSTRT * 0, 

SHOa » 1, 

IRPEAT a 1, 

NZONES a 1, 

ECRAT»3, ASUB»3,2,1.5, NASUB»3, 

ASUP a 5,10,20,30,40,50,60,77.5, NASUPaS, 

RSI a 5.1527, RWTU a 1, rwTO a ,392, THETAI = 25.4167, RI = 1.73921, 
IWALL»4, THE=5.3748, TH£TAa37, 

NWS a 18, 


RS * 0, 

1.. 126405, 

1.173862, 

1 .255959, 

1.379498, 

1.551987, 


1.741611, 

1 .947049, 

2.186550, 

2.212090, 

2.226684, 


2.363046, 

3.785419, 

5.338208, 

6.586136, 

7.564615, 


8.298253, 

8.803421, 




ZS * 0, 

.2990952, 

.3626237, 

.4733999, 

,6422291 , 

.8823416, 


1.152396, 

1.452449, 

1.812502, 

1.851582, 

1.873954, 


2.085214, 

4.546562, 

7.947338, 

11.54691, 

15.34513, 


19.34211, 

23.52318, 




lOFF s 3, 

r 





SEND 






REACTANTS 






H 2. 




100. 

•2154. 1 

0 2. 




100. 

•3102. 1 


ORIGINAL PACE 13 
OF POOR QUALITY 


Thrust Chamber 
Geometry 


20.27F 

90.180 


Tank Conditions 
CPIA 246 


NAMELISTS 
SOOE 
RKT a T, 

P a 3285., XP a 1, PSIA a T, DELH a 0, 
OFSKEO a 6.05485, 

SEND 

REACTIONS 



H * OHa H20 

,M1,Aa7.5E23 ,Na2.6,Ba0., 

(AR) NO. 1 


0 * H a OH 

,M2,Aa4.0E18 ,N>1. ,BaO., 

(AR) NO. 2 


0 ♦ 0 a 02 

,M3,Aa1.2E17 ,Na1, ,8=0., 

(AR) NO. 3 


H * H a H2 

,M4,A«6.4E17 ,N»1. ,B=0., 

(AR) NO. 6 

END TBR REAX 



H2 OH 

a H ♦ H20 

,Aa2.19E13,N=0. ,Ba5.15, 

8AULCH(68)L2 

OH OH 

a 0 H20 

,A«5.75E12,N=0. ,B=.78, 

BAULCH(68)L2 

H OH 

a 0 ♦ H2 

,A=7.33E12,Na0. ,8=7.3, • 

8AULCH(68)L2 

0 ♦ OH 

a H ♦ 02 

,Aa1.3E13 ,N=0. ,8=0., 

BAULCH(69)L3 

LAST .RE AX 




THIRD 300'^ REAX RATE RATIOS 




H1=20'h2o,5'H2,5*O2, 12. 5*h, 12.5*0, 12.5*OH, 
M2a5*M20,5*H2,5*02,12.5'H, 12.5*0,12.5*0H, 
M3a5*H20,5*H2,4.5*02, 12. 5'h, 12.5*0, 12.5*OH, 
M4a20*H2O,4*H2, 1 .5*02,25*H,25*0,25*0H, 

UST CARD 
SODK 

JPRNT*-2, 

EP a 77.5, 

SEND 
STRANS 
HP a 250, 

XM a 1, 

SEND 
SMOC 
NC a 0, 


\ 


ISHCKal, 


Finite Rate 
Che.mistry 


Regen Cooling 
Circuits 


IMA;;a40, IMAXFal, 

SEND 
SBLM 

OFC a 2, 2, DISTRB a 0, 0, 

XCO a -2.4842, 1.85, 

XCE a 1.849, 23.5, 

IHFLAG a 0, 

NTQW a 27, ’ 

XTOW a -2.4842, -2.1348, -1.9407, -1.7467, -.9704, -.7763, 
-.1941, .1941, .3881, .9704, 1.5526, 1.9407, 1.9601, 
5.8222, 6.7926, 7.7629, 8.7333, 9.7037, 11.6444, 13l585. 
15.5258, 17.4666, 19.4073, 21.3480, 23.523178, 

TQW a 1360, 1500, 1510, 1500, 1470, 1470, 1490, 1490, 840, 

820, 790, 760, 1450, 1260, 1060, 960, 890, 850, 830, 

765, 745, 730, 720, 715, 710, 

NTR a 10, 

XINO a -2.4842, -2.40, -2.30, 

RING a 3*1.73251, 

UEO a 535, 690.8, 979.2, 

TEO a 6499.7, 6496.8, 6489.5, 

PEO a 2918.1, 2906.3, 2877, 


3881, 

3.8815, 


830, 

795, 


NSEGS a 6, 

XSEG a -2.5, -.4658, .4658, 
APROF a 50, 76, NPROF « 0, 
KDTPLT»0,KMTPLTaO,KTWPLTaO, 
SEND 


4.6578, 9.3155, 16.3021, 23.53, 


Wall 

Te.mperature 


& 

Figure 3: TDK Input for SSME Current Design Analysis 
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TITLE SSME « LARGE THROAT DESIGN, COOLED WALL 
DATA 
SOATA 

OOE * t, OOK * 1, TO< » 1, BLM * 1, IRSTRT « 0, 

SHOCK a 1, 

IRPEAT a 1, 

NZONES a 1, 

ECRAT=2.7, ASU8*2.7,2,1.5, NASUB=3, 

ASUP a 5,10,20,30,40,50,60,69.5, NASUP=8, 

RSI a 5.4417, RWTU a .494, RWTD a .2, THETAI = 25.4167, RI * 1.73921, 
lUALLa4, 

RS =1.066583 , 1.111520 , 1.189257 

1.306235 , 1.469563 , 1.649117 , 1.843644 , 2.070426 

2.094609 , 2.108428 , 2.237548 , 3.584382 , 5.054705 

6.236357 , 7.162870 , 7.857546 , 8.335885 

2S aO. 1560446 , 0.2161992 , 0.3210922 


0.4809552 
1.626081 
10.80651 
NWS a 

THE a 

THETA a 
lOFF a 3, 
SEND 

REACTANTS 

H 2. 

0 2 . 


0.7083157 
, 1.647265 

, 14.40301 

17, 

5.374800 

37, 


0.9640279 
1 .847305 
18.18771 


1.248146 

4.177936 

22.14673 


100 . 

100 . 


■2154. 

-3102. 


1.589077 

7.398101 


20.27F 

90.180 


NAMELISTS 

SOOE 

RKT a T, 

P a 3010., XP a 1, PSIA a T, DELH 
OFSKED a 6.05485, 

SEND 

REACTIONS 


0 , 


H + OH« H20 


,M1,A=7.5E23 ,N.2.6,B*0., 

(AR) MO. 1 

0 > H a OH 


,M2,A»4.0E18 ,H=1. ,B=0., 

(AR) NO. 2 

0 0 a 02 


,M3,A»1.2E17 ,M.1. ,BaO., 

(AR) NO. 3 

H + H a h2 


,H4,A»6.4E17 ,M«1. ,B=0., 

(AR) NO. 6 

END TBR REAX 




h2 * oh * H ♦ 

H20 

,A»2.19E13,M»0. ,B*5.15, 

BAUICH(68)L2 

OH ♦OH = 0 *► 

H20 

,A*5.75E12,N«0. .B..78, 

BAULCH(68)L2 

H ♦ OH a 0 ♦ 

H2 

,Aa7.33E12,M«0. ,8*7.3, 

BAULCH(68)L2 

0 ♦ OH a H ♦ 

02 

,A*1.3E13 ,M*0. ,B»0., 

BAU(.CH(69)L3 


LAST REAX 

THIRD BODY REAX RATE RATIOS 

Ml =20*H20 , 5*H2 , 5*02 , 1 2 . 5*H , 1 2 . 5*0 , 1 2 . 5*OH , 

M2=5*H20 , 5*H2 , 5*02 , 1 2 . 5*H , 1 2 . 5*0 , 1 2 . 5*0H , 

M3=5*H20,5*H2,4.5*02,12.5*H, 12.5*0, 12. 5*0H, 

M4=20*H20 , 4*H2 , 1 . 5*02, 25*H , 25*0, 25*0H , 

LAST CARO 
SOOK 

JPRNT=-2, 

EP a 69.5, 

SEND 
STRANS 
MP a 250, 

XM a 1, 

SEND 

NCaO, 

1MAX=40, I MAX Fa 1, 

ISHCKa2, 

SEND 

SBLM 

IHFLAG a 0, 

OFC a 2 , 2, OISTRB a 0, 0, 

XCO a -2.4842, 1.85, 

XCE a 1.849, 22,15, 

NTR a 10, 

XINO a -2.5, -2.40, -2.30, 

R I NO a 3*1.6432, 

UEO a 535, 690.8, 979.2, 

TEO a 6499.7, 6496.8 , 6489.5, 

PEO a 2918.1, 2906.3, 2877, 

NSEGS a 6, 

4.6578, 9.3155, 16.3021, 22.16, 

APROF a 50, 76. NPROF a 0, ' ' 

KDTPLTaO,lCMTPLTaO,KTWPLT=0, 

NTQW a 27, 

XTQW a -2.4842, -2.1348, -1.9407, -1.7467, -.9704, -.7763, -.3881 

’-5526, 1.9407, 1.9601, 3.8815, 
, 7.7629 , 8.7333 , 9.7037, 11.^44,' ISISSS, 

T~. 1^-^666. 19.4073 . 21.3480 , 23.523178, 

TOW * 1360, 1500, 1510, 1500, 1470, 1470, 1490, 1490, 840, 830 
820. ^2' ^^50, 1260, 1060 , 960 , 890 , 850 , 830 , 795, 

765, 745, 730, 720, 715, 710, 

SEND 


TDK Input for SSF.E "Large Throat" Design Analysis 
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Figure 4: 
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TDK PERFORMANCE SUMMARY : SSME • COOLED UALL , CURRENT DESIGN, 109X 

REAL UALL CONTOUR 1 ZONES UALL TEMP 


FIRST 

TOK/BLM SOLUTION 


SECOND 

TDK/BLM SOLUTION 


THRUST CHAMBER OPERATING CONDITIONS 


CHAMBER PRESS 

CPSIAl 

3285.000 

CHAMBER TEMP 

CR] 

6543,192 

MIXTURE RATIO 

(•] 

6.054850 

H (OXID) 

(CAL/MOLE3 

*3102.000 

H (FUEL) 

(CAL/MOLEl 

*2154,000 

HCHAM (ODE) 

CBTU/LB] 

•422.3787 

DELH (AVERAGE) 

[BTU/LB1 

O.OOOOOOOE^OO 

DELHI (AVE) 

CBTU/LB] 

O.OOOOOOOE^OO 

THRUST CHAMBER GEOMETRY 


ECRAT 

C-3 

3.000000 

RI 

C-3 

1.739210 

THETAI 

(DEGREES! 

25.41670 

RUTU 

C-3 

1.000000 

RSTAR 

C1NCHES3 

5.152700 

RWTO 

C-3 

0.3920000 

NIT 

C-3 

243.0000 

THE 

(DEGREES! 

5.374800 

THETA 

(DEGREES! 

37.00000 

EP (NOZZLE) 

C-3 

77.50022 

SHOCK PARAMETERS 


XA (X/R*) 

C-3 

0.2359115 

XB (X/R*) 

C-3 

O.OOOOOOOE^O 

SHOCK (R/R*) 

C-3 

2.584894 

SHOCK (X/R*) 

C-3 

8.834632 

SHOCK STRENGTH 

C-3 ^ 

0.2892100 


EXIT FLOU PROPERTIES 


P (AXIS, EXIT) 

CPSIAl 

0.2726767 

P (WALL, EXIT) 

CPS I A] 

6.109700 

T (WALL,EXIT) 

CR] 

2483.953 

V (WALL, EXIT) 

CFT/SEC3 

14039.21 

HA (WALL,EXIT) 

1 C-3 

4.271464 


ONE-OI^ENSIONAL FLOU PERFORMANCE 

ISP (ODE) 

(SECONDS! 

465.1311 

ISP (OOK) 

(SECONDS! 

464.9658 

ISP (OOF) 

[SECONDS! 

446.7661 


TUO-OIHENSIONAL FLOW PERFORMANCE 

CO 

C-3 

0-9896466 

CF (TDK) 

C-3 

1.921868 

CSTAR (TDK) 

CFT/SEC3 

7713.456 

THRUST (TDK) 

(POUNDS! 

526597.2 

WOOT (TDK) 

CLB/SECl 

1142.907 

ISP (TDK) 

(SECONDS! 

460.7523 


BOUNDARY LAYER PARAMETERS 


OFOPT (8LME) 

[POUNDS] 

8386.935 

DF (8LME) 

(POUNDS! 

7931.962 

OISP (BLHE) 

(SECONDS! 

6.940162 

THETA (EXIT) 

(FEET! 

0.1777533E*01 

DEL* (EXIT) 

(FEET! 

0.2187267E-01 

EP (REGEN) 

(BTU/SEC! 

77.46112 

SCOOT (REGEN) 

(BTU/SEC! 

135682.6 

SQOOT (LOSS) 

(BTU/SEC! 

56.39063 


THRUST CHAMBER PERFORMANCE 


THRUST (TO 

(POUNDS! 

518665.2 

ISP (TO 

(SECONDS! 

453.3121 

CF (TO 

(-3 

1.892919 


3285. 000 
6609.458 
6.054850 
•3102.000 
-2154.000 
•303.6657 
O.OOOOOOOE-fOO 
118.7092 


2.997371 

1,739210 

25.41670 

0.9985945 

5.156323 

0.3910218 

243.0000 
5.299142 

37.00000 
76.50255 


0.2359115 

O.OOOOOOOE+00 

2.589131 

8.839259 

O.OOOOOOOE4-00 


0.2756312 

6.161937 

2561.766 

14149.57 

4.243621 


469,2740 

469.0910 

449.4862 


0.9910184 
1.928235 
7768.724 
528341 ,6 
1136.373 
464.9368 


8285.425 

7829.081 

6.889535 

0.1777533E-01 

0.2187267E-01 

77.46112 

O.OOOOOOOE-^00 

56.39063 


520512.5 

458.0472 

1.899662 


Figure 6 


TDK 
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TDK PERfORMAMCE SUMMARY : SSME • LARGE THROAT DESIGN, COOLED WALL 

REAL WALL CONTOUR 1 ZONES WALL TEMP 
INDUCED SHOCK AFTER XA 


FIRST 

TOK/BLM SOLUTION 


THRUST CHAMBER OPERATING CONDITIONS 


CHAMBER PRESS 

CPS I A] 

3010.000 

CHAMBER TEMP 

CR] 

6526.882 

MIXTURE RATIO 

C-] 

6.054850 

H COXIO) 

CCAL/MOLE) 

-3102.000 

H (FUEL) 

CCAL/MOLE] 

-2154.000 

HCHAM (ODE) 

CBTU/LB] 

-422.3749 

DELH (AVERAGE) 

CBTU/LB] 

O.OOOOOOOE+00 

DELHI (AVE) 

CBTU/IB] 

THRUST CHAMBER GEOMETRY 

O.OOOOOOOE+00 

ECRAT 

C-] 

2.700000 

R1 

C-] 

1.739210 

THETA I 

[DEGREES) 

25.41670 

RWTU 

C-1 

0.4940000 

RSTAR 

CINCHES] 

5.441700 

RUTO 

C-] 

0.2000000 

NIT 

C-] 

244.0000 

THE 

[DEGREES] 

5.374800 

THETA 

[DEGREES] 

37.00000 

EP (NOZZLE) 

[-] 

SHOCK PARAMETERS 

69.48698 

XA (X/R*) 

[-1 

0.1203630 

XB (X/R*) 

C-] 

O.OOOOOOOE-^00 

SHOCK (R/R*) 

C-1 

2.637427 

SHOCK (X/R*) 

(■] 

8.547814 

SHOCK STRENGTH 

1 C-] 

EXIT FLOW PROPERTIES 

1.131168 

P (AXIS, EXIT) 

CPSIA] 

0.2820055 

P (WALL, EXIT) 

CPSlAl 

6.195770 

T (WALL, EXIT) 

CR] 

2519.777 

V (WALL, EXIT) 

EFT/SEC] 

13992.76 

MA (WALL, EXIT! 

1 C-3 

ONE-DIMENSIONAL- FLOW PERFORMANCE 

4.229021 

ISP (ODE) 

[SECONDS] 

463.2604 

ISP (OOK) 

CSECONOS] 

463.0793 

ISP (OOF) 

[SECONDS] 

TWO-DIMENSIONAL FLOW PERFORMANCE 

444.7228 

CO 

C-] 

0.9802917 

CF (TDK) 

[-] 

1 .892262 

CSTAR (TDK) 

CFT/SEC] 

7774.756 

THRUST (TDK) 

(POUNDS] 

529866.9 

WOOT (TDK) 

[LB/SEC] 

1158.788 

ISP (TDK) 

[SECONDS] 

BOUNDARY LAYER PARAMETERS 

457.2597 

DFOPT (8LME) 

[POUNDS] 

8341.068 

OF (BLME) 

[POUNDS] 

7889.169 

DISP (BLME) 

[SECONDS] 

6.308123 

THETA (EXIT) 

[FEET] 

0.1780233E-01 

DEL* (EXIT) 

[FEET] 

0.2142398E-01 

EP (REGENT 

CBTU/SECl 

69.48714 

SQDOT (RE.EN) 

[8TU/SEC] 

135283.0 

SQOOT (LOSS) 

[BTU/SEC] 

THRUST CHAMBER PERFORMANCE 

o,oooooooE-*-oo 

THRUST (TO 

[POUNDS] 

521977.8 

ISP (TO 

[SECONDS] 

450.4516 

CF (TO 

[-1 

1.864088 


Figure 7: TDK Performance Summary for SSME 

Configuration 


SECOND 

TDK/BLM SOLUTION 


3010.000 

6591.104 

6,054850 

-3102.000 

-2154.000 

-305.6916 

O.OOOOOOOE+00 

116.6856 


2.694912 

1.739210 

25.41670 

0.4925624 

5.446941 

0.1988453 

244.0000 
5.327114 

37.00000 
68.57317 


0.1203630 

O.OOOOOOOE+00 

2.612652 

8.626249 

1.138703 


0.2786561 

6.247626 

2589.486 

14109.94 

4.210537 


467.2708 

467.0724 

447.3469 


0.9821196 
1 .898804 
7830.471 
531698.8 
1152.760 
461.2397 


8277.303 

7824,075 

6.787252 

0.1780233E-01 

0.2142398E-01 

69.48714 

O.OOOOOOOE+00 

O.OOOOOOOE+OO 


523874.7 
454.4525 
1 .870862 
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21 


>(ai MO 

r IV O'ii H i B 1 ^VJ s s 


LJl 

cn 


ORIGINAL PAGE IS 
OF POOR QUALITY 




( • 

1 # # / 1 

t » » 

» 1 • » 

t 1 1 1 « 1 

r * » 

o 

S3 ^ 

s ? 

o 

— C; *i 

O 3* 

c 

o 

o 


o 

o 

o 

o 


V 

Cjbr 

Oo 





CO 




:r 

I!2 


rO 




I 

o 

o 


Ln 



> 


22 
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SENSITIVITY OF HASS FLOW RATE TO FLOW ANGLE 
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Plj-ure 13: Sensitivity of Nozzle Mass Plow Coefficient, C^, to Streamline 

Angle Variation 



COnPARISON OF TDK (MOD. SAUER) WITH VNAP2 
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Purdue Results 
SEA Inc. Results 



Figure 15: SSME Current Design, Induced Shock Predictions 







Large Throat” Design, Induced Shock Prediction 













starting Conditjlonsi 



Figure 19: Shock Strength vs x/r for Differing Transonic Start Lines 
"Large Throat" SSME Design. 






SUMMARY OF SHOCK STRENGTH PREDICTIONS 
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X/R* = 9 IS SHOCK POINT POSITION ON LAST LRC. 
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APPENDIX A 


DISCUSSION OF SUMMARY OUTPUT 

At the end of each computer run, the TDK program prints a 
table of summary output. Examples of this output are given in 
Figures 6 and 7. The first item printed is the title 

TDK PERFORMANCE SUMMARY: (title) 

Next, performance parameters summarizing the results of the 
calculations are printed in three columns. The left hand column 
identifies each item to be printed and its units. The first 
column of results is for the first MOC/BLM solution, and second 
column of results is for the second MOC/BLM solution. 

The items labeled ISP(TDK), THRUST(TDK), WDOT(TDK), CD, 
CSTAR(TDK), and CF(TDK)' are defined as: 


N ame 

Units 

Definition 


ISP(TDK) 

sec 

Specific Impulse, I 

sp , 2D 

THRUST(TDK) 

Ibf 

Thrust, Fgjj 


WD0T(TDK) 

Ibm/ sec 

Mass Flow, thgu 


CD 

- 

Flow coefficient, C 

D , 2D 

CSTAR(TDK) 

ft/ sec 

Characteristic Velocity, 

CF(TDK) 

- 

Thrust coefficient. 

^F,2D 

If a TDF or 
these items will 

TDE calculation 
be relabeled. 

was made rather than 

TDK, then 
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The values given for and are obtained by 
evaluating the integrals discussed in the documentation for MOC 
module subroutine CHAR. The mass flow integral is evaluated 
across the MOC initial data line. The integral is evaluated 

3 p 

across the initial data line, plus the integral of the axial 
component of wall pressure from the initial data line to the 
nozzle end point. The thrust is 


2D 


^sp,2D ”2D 


The nozzle flow coefficient is calculated as 

C. 

where 


'D,2D 


*2D^ *1D 


A 


1 D 


« « » 

PiD ^D A 


is obtained from ODE. A is the geometric throat area 

« 

The characteristic velocity, C , is calculated as 

* 


'2D 


^Sp,2D ®'^*^F,2D 


where the thrust 

coefficient 

, Cp, is defined 

by 


^F,2D 

- F^„ / P A* 
2D c 


» 

Again, A , is the geometric 

throat area. 

• 

The items 

labeled DISP(BLME), DFOPT(BLME), DF(BLME) 

ISP (TO, THRUST( 

TC), and CF(TC) are defined 

as follows: 

N ame 

Units 

Definition 

• 

DISP(BLME) 

sec 

^^sp,BLM 


DFOPT(BLME) 

Ibf 


^amb 

DF( BLME) 

Ibf 

^^BLM 
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Name 


Units 

Definition 

ISP(TC) 


sec 

^sp,2D,BLM 

THRUST(TC) 

Ibf 

^2D,BLM 


CF(TC) 


— 

2D, BLM 


In 

the above. 

the thrust loss calculated 

by the 

following 

relation: 



A F * 2 ir r 

■ ^ BLM e 

cos * 

e 

Pe ^e ®BLM 



- 2irr 

e 

cos a * 
e 

- %Mb> « 

* 

BLM 

The additional items, printed are 


^^sp,BLM 

m 

^^BLM^ "*2D 


^sp 

, 2D, BLM. 

m 

(F2d“ ^^bLM^ 


'.2D 

,BLM 

m 

^2D * ^^BLM 


= F. 

2D, BLM 

m 

^2D,BLM ^ ^c 

» 

A 


the BLM is given 


Per formance. parameters summarizing the results of the second 
MOC/BLM calculation differ from the first set printed as 
follows : 

1) The TDK and ODE results for the second (invisid 

core) run have a modified wall contour, throat area, 
and expansion ratio. Note that the system enthalpy is 
different than used in the first run if the BLM wall 
calculation is not adiabatic. 

2) The BLM results are re-calculated using edge conditions 
from the second TDK run. BLM itself is not re-run, and 

* 

the values of 6 and 9 are not changed. 

The equations used are the same as given above. 
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